home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2001 May / SGI IRIX Base Documentation 2001 May.iso / usr / share / catman / p_man / cat3 / complib / ssytri.z / ssytri
Encoding:
Text File  |  1998-10-30  |  2.8 KB  |  133 lines
  1.  
  2.  
  3.  
  4. SSSSSSSSYYYYTTTTRRRRIIII((((3333FFFF))))                                                          SSSSSSSSYYYYTTTTRRRRIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SSYTRI - compute the inverse of a real symmetric indefinite matrix A
  10.      using the factorization A = U*D*U**T or A = L*D*L**T computed by SSYTRF
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE SSYTRI( UPLO, N, A, LDA, IPIV, WORK, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, N
  18.  
  19.          INTEGER        IPIV( * )
  20.  
  21.          REAL           A( LDA, * ), WORK( * )
  22.  
  23. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  24.      SSYTRI computes the inverse of a real symmetric indefinite matrix A using
  25.      the factorization A = U*D*U**T or A = L*D*L**T computed by SSYTRF.
  26.  
  27.  
  28. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  29.      UPLO    (input) CHARACTER*1
  30.              Specifies whether the details of the factorization are stored as
  31.              an upper or lower triangular matrix.  = 'U':  Upper triangular,
  32.              form is A = U*D*U**T;
  33.              = 'L':  Lower triangular, form is A = L*D*L**T.
  34.  
  35.      N       (input) INTEGER
  36.              The order of the matrix A.  N >= 0.
  37.  
  38.      A       (input/output) REAL array, dimension (LDA,N)
  39.              On entry, the block diagonal matrix D and the multipliers used to
  40.              obtain the factor U or L as computed by SSYTRF.
  41.  
  42.              On exit, if INFO = 0, the (symmetric) inverse of the original
  43.              matrix.  If UPLO = 'U', the upper triangular part of the inverse
  44.              is formed and the part of A below the diagonal is not referenced;
  45.              if UPLO = 'L' the lower triangular part of the inverse is formed
  46.              and the part of A above the diagonal is not referenced.
  47.  
  48.      LDA     (input) INTEGER
  49.              The leading dimension of the array A.  LDA >= max(1,N).
  50.  
  51.      IPIV    (input) INTEGER array, dimension (N)
  52.              Details of the interchanges and the block structure of D as
  53.              determined by SSYTRF.
  54.  
  55.      WORK    (workspace) REAL array, dimension (N)
  56.  
  57.      INFO    (output) INTEGER
  58.              = 0: successful exit
  59.              < 0: if INFO = -i, the i-th argument had an illegal value
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. SSSSSSSSYYYYTTTTRRRRIIII((((3333FFFF))))                                                          SSSSSSSSYYYYTTTTRRRRIIII((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              > 0: if INFO = i, D(i,i) = 0; the matrix is singular and its
  75.              inverse could not be computed.
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.